A Redshift Determination for XRF 020903: First Spectroscopic 

Observations of an X-Ray Flash 



A. M. Soderberg, S. R. Kulkarni, E. Berger, D. B. Fox, P. A. Price, S. A. Yost and M. P. 

Hunt 

Division of Physics, Mathematics and Astronomy, MS 105-24, California Institute of 

O ' Technology, Pasadena, CA 91125 

O ' 

^ : D. A. Frail and R. C. Walker 

O 

Q ■ National Radio Astronomy Observatory, P. 0. Box 0, Socorro, NM 87801 
^ ' M. Hamuy and S. A. Shectman 

m 

^ \ Carnegie Observatories, 813 Santa Barbara Street, Pasadena, CA 91101 

O . and 



o 
o 



X 



J. p. Halpern and N. Mirabal 

Astronomy Department, Mailcode 5246, Columbia University, New York, NY 10027 

ABSTRACT 

We report the discovery of optical and radio afterglow emission from the 
extremely soft X-ray flash, XRF 020903. Our spectroscopic observations provide 



H ' the first redshift for an X-ray flash, thereby setting the distance scale for these 

events. At z = 0.251, XRF 020903 is one of the nearest cosmic explosions 
ever detected, second only to the recent GRB 030329 and the unusual GRB 
980425/SN 1998bw. Moreover, XRF 020903 is the first X-ray flash for which 
we detect an optical afterglow. The luminosity of the radio afterglow of XRF 
020903 is 1000 times greater than that of Ibc supernovae but similar to those of 
GRB afterglows. From broadband afterglow modeling we show that the explosion 
energy of XRF 020903 is not dissimilar from values inferred for typical gamma-ray 
bursts, suggesting that these cosmological explosions may derive from a similar 
mechanism. 



Subject headings: gamma-ray bursts: specific (XRF 020903) — supernovae: gen- 
eral 
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1. Introduction 

Prior to the detection of afterglows, gamma-ray bursts (GRBs) were enshrouded in mys- 
tery for nearly thirty years. Great progress in our understanding of these energetic events 
came with the first redshift measurement which placed GRBs at cosmological distances (Met- 
zger et al. 1997). In a similar fashion, the mystery of X-ray flashes (XRFs) has been fueled 
by the absence of confirmed redshifts. These events, identified in the 1990s by BeppoSAX, 
are characterized by a peak energy in i/F^, of £^peak ~ 25 keV (compared to Epeak ~ 250 
keV for GRBs). With a distribution of durations similar to those observed for GRBs, it has 
been assumed that XRFs are associated with GRBs and therefore share their extragalactic 
distance scale (Heise et al. 2001). 

The subsequent discovery of XRF X-ray and radio afterglows with properties similar to 

those observed in GRB afterglows has further strengthened this association (Harrison et al. 
2001; Amati et al. 2002; Taylor, Frail & Kulkarni 2001). Still, the question of whether the 
difference between GRBs and XRFs is intrinsic or extrinsic remains unanswered. Assuming 
XRFs are simply GRBs observed away from the jet collimation axis, they would have less 
7-ray emission and the difference is extrinsic, based solely on the line-of-sight to the observer. 
On the other hand, the difference could be intrinsic, namely XRFs may represent a class 
of explosions which are similar in energetics to GRBs yet characterized by less relativistic 
ejecta possibly due to a heavier baryonic load. It is clear that by setting the distance scale 
for XRFs (and hence their energy scale) we can begin to distinguish between extrinsic and 
intrinsic effects. 

In this paper we present the first spectroscopic redshift for an X-ray flash, XRF 020903, 
placing this event among the nearest high energy explosions and offering confirmation that 
X-ray flashes are cosmological and produce a total energy output similar to that observed 
in GRBs. 



2. Observations 

On 2002 September 3.421 UT the Wide-Field X-ray Monitor (WXM) and Soft X-ray 
Camera (SXC) aboard the High Energy Transient Explorer-2 (HETE-2) detected an X-ray 
flash within the 0.5-10 keV energy band. With an exceptionally low peak energy of -Epcak ~ 5 
keV and a fluence of 7.2 x 10^^ erg cm~^, XRF 020903 is the softest event ever detected 
by HETE-2 with a ratio of X-ray fluence (Sx) to 7-ray fluence (S^) of \og{Sx/ S^^) = 4.3 
(Sakamoto et al. 2003). Ground analysis provided a localization for XRF 020903 centered at 
a(J2000) = 22''49°'0P, (5(J2000) = -20°55'47" with a 4' x 31' uncertainty region at At ^ 0.3 
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days (Ricker et al. 2002). 

2.1. Ground-based Photometry 

We began observing the field of XRF 020903 on 2002 September 4.32 UT {At 0.9 
days) with the Palomar Observatory 200-inch telescope (P200) equipped with the Large Field 
Camera (LFC). With a total exposure time of 20 minutes under photometric conditions 
(stellar FWHM ~ 1.2") the observations reached a limiting magnitude of i? ~ 23 mag. 
Visual comparison with Digitized Sky Survey archival images did not reveal an afterglow 
candidate. 

A second epoch was obtained on September 10.30 UT {At Ri7 days) using the same 
observational set-up and in similar observing conditions. Image subtraction between the first 

and second epochs revealed one variable object within the HETE-2 error region (Figure 1) 
located at a(J2000) = 22'H8°'42.34% 5(J2000) - 20°46'09.3" and lying ~ 4 arcsec NW of a 
bright elliptical galaxy (hereafter G2; Soderberg et al. 2002). 

With an approximate magnitude of i? ~ 19 at At ?s 1 day, the new object decreased 
in brightness by 1.4 magnitudes between the two epochs, implying a temporal flux decay 
index of a ~ — 1. As the source proved consistent with a typical GRB afterglow evolution, 
the optical transient was adopted as a suitable candidate for the optical afterglow of XRF 
020903. 

We observed the afterglow position on three additional epochs with the P200 and the 
MDM Observatory 1.3 meter telescope (Table 1 and Figure 2). These, along with Digitized 
Sky Survey archival images, reveal the presence of an extended source with i? ~ 21 mag 
coincident with the position of the transient and thereby suggestive of a host galaxy (hereafter 
Gl). 

Due to the underlying extended source and the proximity of the transient to G2, accurate 
magnitude estimates relied on careful PSF photometry of the field. Absolute calibration of 
field stars was supplied by Henden (2002). We used 12 unsaturated field stars in common 
between the Palomar and Henden images to fully calibrate observations of the variable source. 
We note that although the transient decreased in brightness quickly between the first two 
epochs, later epochs {At 30 — 40 days) appear to indicate a plateau in the fight curve 
which was confirmed by other observers (Gorosabel et al. 2002) and may originate from 
unresolved fiux contamination from Gl. 
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2.2. Ground- based Spectroscopy 

Initial spectroscopy of the transient was performed on 2002 September 28.1 UT with 
the Magellan Baade Telescope using the Low Dispersion Survey Spectrograph (LDSS2). A 
position angle of 168 degrees was used such that spectral information on Gl and G2 was 
obtained simultaneously. Despite the positional coincidence of G2 and Gl, it was found 
that the systems are not physically associated, separated by 3900 km/s in velocity space 
(Soderberg et al. 2002). In this epoch, the transient source was still significantly bright. 

Further spectroscopic observations were made of the putative host (Gl) with the Echelle 
Spectrograph & Imager (ESI) mounted on the Keck II telescope on 2003 July 4 UT (At f« 300 
days). These observations do not include any flux from the transient source. During a total 
exposure time of 1 hour, we obtained a spectrum of Gl and G2 with a slit width of 0.75 
arcsec. 

We flnd that G2 is a large eUiptical galaxy with at least one interacting galaxy com- 
panion, G3, located < 1" to the NE (see Figure 4). Observations of G2 exhibit a relatively 
smooth continuum with features typical of an elliptical galaxy. The Ca II H and K absorp- 
tion lines give a redshift oi z — 0.235, while G3 is offset by only 240 km/s at a redshift of 
z = 0.236. 

In contrast, Gl is shown to be an active star-forming galaxy at z=0.251 with a rich set of 
narrow bright emission lines (Figure 2). The [OIII]/H/? and [Nil] /Ha intensity ratios indicate 
that the galaxy is a low-metallicity and high-excitation starburst galaxy. In addition, the 
flux ratio of the [Nelll] and [Oil] lines is ^^3869^^3727 0.43, similar to the observed value 
for the host galaxy of GRB 970508 (Bloom et al. 1998) and approximately 10 times higher 
than typical values for H II regions. The bright [Nelll] emission lines observed in GRB hosts 
are thought to be indicative of a substantial population of massive stars. On the other hand, 
Chornock & Filippenko (2002) note that a spectrum of Gl taken with the Low Resolution 
Imaging Spectrometer (LRIS) on the Keck I telescope reveals a deflcit of emission at rest 
wavelengths < 4000A which is consistent with a population of older stars. 

2.3. Hubble Space Telescope 

The afterglow candidate was observed with the Hubble Space Telescope (HST) using 
the Advanced Camera for Surveys (ACS) under Program No. 9405 (RL: A. Fruchter). Three 
epochs of imaging were obtained from At ~ 94 to ~ 300 days with exposure times of 1840 
sec in the F606W fllter. Following "On-The-Fly" pre-processing the data were drizzled using 
standard IRAF tools (STSDAS; Fruchter & Hook 2002). In drizzhng our flnal images, we 
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retained the native WFC pixel scale of 0.05" and used a pixf rac of 1.0. The HST images 
reveal a complex galaxy morphology for Gl, suggesting a system of at least four interacting 
galaxies (see Figure 4). 

To locate the optical transient with respect to the host galaxy complex, we performed a 
source to source comparison of our first-epoch Palomar (2002 September 4) and HST (2002 
December 3) images. We found 42 unsaturated, unconfused sources in common between 
these two images, and were able to match the two coordinate lists with an rms mapping 
uncertainty of 0.06 arcsec. We derived the position of the transient from the difference 
image of the 2002 September 4 and September 10 data. Since the optical transient is well- 
detected in the difference image, the uncertainty in its centroid position is negligible. The 
uncertainty in the coordinate mapping thus dominates the uncertainty in the position of the 
source relative to the host galaxy complex. The optical transient appears to overlap with 
the SW component of the Gl complex (Figure 4). 

2.4. Radio Observations of XRF 020903 

2.4-1- Very Large Array Data 

We began observations of the field of XRF 020903 with the Very Large Array (VLA^) 
on 2002 September 27.22 UT. A radio source was detected in coincidence with the optical 
transient at a location of a(J2000) = 22M8°^42.34% 5( J2000) = -20°46'08.9" with an uncer- 
tainty of 0.1 arcsec in each coordinate. The initial observation showed the radio source to 
have a flux density of Fiy = 1.06 ± 0.02 mJy at 8.5 GHz. The National Radio Astronomy 
Observatory VLA Sky Survey (NVSS; Condon et al. 1998) did not show any evidence for 
a pre-existing source at this location down to a limit of 1 mJy. Further observations at 8.5 
GHz on September 29.11 {At «i 26 days) showed that the source faded to Fj, = 0.75 ± 0.04 
mJy. We continued monitoring the transient source with the VLA over the next ^ 370 days 
at frequencies of 1.5, 4.9, 8.5 and 22.5 GHz. The lightcurve is displayed in Figure 5. 

2-4-2. Very Long Baseline Array Data 

The relatively low redshift and strong radio emission of the transient source made it an 
ideal candidate for Very Long Baseline Array (VLBA) observations. On September 30.23 



^The VLA is operated by the National Radio Astronomy Observatory, a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, Inc. 
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UT {At 27 days) we observed the radio transient for a total duration of 6 hours at 8.5 
GHz. Using sources J2253+1608 and J2148+0657 wc were able to flux calibrate the field 
and we phase referenced against source J2256-2011 at a distance of < 2 degrees. Data 
were reduced and processed using standard VLBA packages within the Astronomical Image 
Processing System (AIPS). We detected the radio transient with a flux of 0.89 ± 0.17 mJy 
at a position coincident with the VLA and optical observations at a location of a(J2000) = 
22h4gm42 33912^ ± 0.00003, (5(J2000) = -20°46'08.945" ± 0.0005. This is our most accurate 
position measurement for the transient object. The source is unresolved within our VLBA 
beam size of 1.93 x 1.73 mas. Furthermore, the consistency between the flux measured with 
VLBA and low resolution VLA observations rules out the presence of diffuse components 
(e.g. jets). 



3. The afterglow of XRF 020903 

The large error box and the proximity of the optical transient to G2 and G3 delayed 
rapid identiflcation of the transient (Price, Schmidt & Axelrod 2002; Pavlenko, Rumyantsev 
& Pozanenko 2002; Uemura et al. 2002; Pruchter et al. 2002). As a result there was little 
optical followup of the transient at early times (see Table 1). At later times, the transient is 
rapidly dominated by the emission from Gl. An additional complication was introduced by 
Gal- Yam (2002) who, based on archival photographic plates from 1954 and 1977, proposed 
that Gl hosted a variable active galactic nucleus (AGN). 

In contrast to the above discussion, the extensive radio light curve (Figure 5) provides 
strong evidence that the transient is the afterglow of XRF 020903. The radio light curve is 
similar to the afterglow of GRBs (Frail et al. 2003). At 4.9 GHz, the source was observed 
to rise to a peak flux at At ~ 29 days and subsequently decay with a characteristic index 
of a —1.1. A steeper dechne of F^, oc t'^-^ was observed at 8.5 GHz where the peak 
flux precedes the flrst observational epoch. We note that the dynamic range of the radio 
afterglow emission clearly distinguishes the source from a radio bright AGN. 

The afterglow interpretation is consistent with the optical spectroscopic observations of 
Gl, namely the absence of broad lines or features typical of AGN, as well as the absence 
of any jet structure in the VLBA images of the radio transient. We proceed accepting the 
notion that the optical/radio transient is the afterglow of XRF 020903. 
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4. Energetics 

At a redshift of z = 0.251, the X-ray flucncc implies an isotropic equivalent energy 
Eiso ~ 1.1 X 10^^ erg. This value is 3 to 6 orders of magnitude lower than the isotropic 
energies of gamma-ray bursts (Frail et al. 2001). Since there is no evidence for a jet break in 
the afterglow observations, the data are consistent with a wide jet or a spherical explosion. 
This, along with the early dechne of the optical flux [F^, oc t~^'^ at At 1 d; Figure 2) rules 
out the possibility that the inferred low energy is due to a viewing angle significantly away 
from the jet axis. In such a scenario, the optical flux would rise until the edge of the jet 
enters the observer's line of sight, peaking at > tjeu where tj^t is the observed time of 
the jet break (e.g. Granot et al. 2002). 

Figure 6 shows that although the energy in the prompt emission of XRF 020903 is 
significantly lower than that observed for typical GRBs, the peak luminosity of the radio 
afterglow is comparable to that found in GRB afterglows. This indicates that XRF 020903 
has a similar kinetic energy to GRB afterglows. To study this in more detail we used standard 
broadband afterglow models with spherical and collimatcd cjccta expanding into circumburst 
media with uniform density and wind density profile, p oc r^^ (e.g. Berger et al. 2000). 
We find that independent of the assumed model, the total kinetic energy is £^ ~ 4 x 10^° 
erg with reasonable values for the circumburst density and energy fractions (electron and 
magnetic field) of n «i 100 cm~^, e^. ^ 0.6 and ^ 0.01, respectively. Thus we confirm the 
afterglow energy is similar to values found for GRB afterglows (Panaitescu & Kumar 2002; 
Berger, Kulkarni & Frail 2003). 

Recent evidence suggests a standard total energy yield (10^^ erg) for all GRBs (Berger 
et al. 2003b), where the total energy yield (-Btot) is defined as the sum of the energy in 
the prompt emission {E.y) plus the mildly relativistic energy as inferred from the afterglow 
(ii^rad) • Clearly, for XRF 020903 the explosion energy is dominated by the mildly relativistic 
afterglow, while a minor fraction of energy couples to the high Lorentz factors characterizing 
the prompt emission. 

Figure 6 also compares the radio luminosity of XRF 020903 with other spherical explo- 
sions - type Ibc supernovae. With a radio luminosity ~ 1000 times greater than typical Ibc 
SNe, it is clear that XRF 020903 is a significantly more energetic explosion. This suggests 
that despite their similar explosion geometries, XRF 020903 and type Ibc supernovae are 
intrinsically different - possibly due to the presense of a central engine (Berger et al. 2003a). 
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5. Conclusions 

We present radio and optical observations of the afterglow of XRF 020903 - the first 
X-ray fiash to have a detected optical afterglow and a spectroscopic distance determination. 
At a redshift oi z = 0.251, this burst has set the distance scale for XRFs and confirmed the 
assumption that they are cosmological in origin. The host galaxy of XRF 020903 appears 
to be a typical star-forming galaxy similar to those of GRB host galaxies. 

The isotropic energy release of the prompt (X-ray) emission of XRF 020903 is at least 
two orders of magnitude smaller than those of GRBs. However, from our broadband modeling 
of the afterglow we determine that the total kinetic energy is about 10^° erg, not dissimilar 
to those inferred for GRBs (Panaitescu & Kumar 2002; Berger, Kulkarni & Frail 2003). In 
comparison with a larger sample of cosmic explosions, XRF 020903 is a clear example where 
less energy is coupled to high Lorentz factors. This source highlights the diversity in high 
energy transients and underscores the importance of studying spherical explosions. 

We thank J. Ulvestad for his help with scheduling the VLBA observation. MH ac- 
knowldges support for this work provided by NASA through Hubble Fellowship grant HST- 
HF-01139.01-A awarded by the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS 
5-26555. JPH and NM acknowledge support from NSF grant AST-0206051. AMS is sup- 
ported by an NSF graduate research fellowship. GRB research at Caltech is supported by 
NASA and NSF. 
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2002 Sep 4.32 UT 2002 Sep 10.30 UT Residual 



Fig. 1. — The field of XRF 020903 was observed with the Palomar 200 inch telescope 
equipped with the Large Field Camera on 2002 September 4.32 and September 10.30 UT. 
Image subtraction techniques revealed one variable object within the HETE-2 error-box 
lying ~ 4 arcsec NW of a bright elliptical galaxy (G2). The residual image clearly indicates 
a transient stellar source which decreased in brightness by 1.4 magnitudes between the two 
epochs. 
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Time Since the Burst (days) 



Fig. 2. — Optical (R band) lightcurve for the transient source. Data were obtained from 
the Palomar 200 inch and MDM 1.3 meter telescopes. The solid line is our best fit to the 
broadband data, including the afterglow component with a temporal decay of a ~ —1.3 
(dashed line) and the host galaxy (Gl) contribution of i? ~ 20.8 mag. The early decline 
rules out the possibility that the viewing angle is significantly away from the jet axis (see 
Granot et al. 2002). 
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Fig. 3. — On 2003 July 4 UT we obtained spectroscopic observations of the putative host 
galaxy (Gl) underlying the optical transient source associated with XRF 020903. Data were 
taken with the Echcllc Spectroscopic Imager (ESI) mounted on the Keck II telescope. The 
source is shown to be an active star-forming galaxy at z=0.251 with a rich set of narrow 
bright emission lines. 
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Fig. 4. — The transient discovered within the error-box of XRF 020903 was observed with 
the Hubble Space Telescope (HST) using the Advanced Camera for Surveys (ACS) on 2002 
December 3 UT. The HST image reveals a complicated galaxy morphology for Gl, suggesting 
a system of at least four interacting galaxies. The location of the optical transient is noted 
on the 7 x 7" HST cutout with a circle which represents a 2a positional uncertainty of 0.12" . 
Nearby galaxies G2 and G3 are labeled accordingly on the HST image. The optical transient 
appears to overlap with the SW component of the Gl complex. 
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Fig. 5. — We observed the field of XRF 020903 with the Very Large Array over the period 
At ^ 25 — 370 days. A bright radio counterpart was detected at a position coincident with 
the optical transient. We continued monitoring the source with the VLA over the next ^ 300 
days at frequencies of 1.5, 4.9, 8.5 and 22.5 GHz. Assuming the transient source is the radio 
afterglow component associated with XRF 020903, we overplot the best fit broadband model 
(solid line) which predicts a total kinetic energy of 4 x 10^° erg in the afterglow. Allowing 
for a host galaxy emission component, we find a better fit assuming a host galaxy fiux of 
-f'i.4GHz ~ 80/iJy (dashed line). The implied star-formation rate from the host is 5 M0yr~^ 
following the conversion from Yun & Carilli (2002). The 8.5 GHz data point marked by 
the open diamond symbol denotes the VLBA observation of the unresolved transient source 
(beam size=1.93 x 1.73 mas). 
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Fig. 6. — The histogram above compares the peak radio luminosity for various cosmic ex- 
plosions, including type Ib/c supcrnovae (spherical ejecta geometry) as well as typical long- 
duration GRBs (coUimated ejecta geometry). Although the energy in the prompt emission 
of XRF 020903 is 10^ - 10^ times lower than that observed for typical GRBs, the total 
kinetic energy in the afterglow of XRF 020903 is comparable to that found in typical GRBs 
afterglows. In comparison with a larger sample of cosmic explosions, XRF 020903 is a clear 
example where there is less energy coupled to high Lorentz factors. 
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Table 1. Optical Observations of the afterglow of XRF 020903. 



Date 


At 


Telescope 


R 


(UT) 


(days) 




(mag) 


2002 Sep 4.32 


0.9 


Palomar 200" 


19.23±0.10 


2002 Sep 10.30 


6.9 


Palomar 200" 


20.60±0.10 


2002 Sep 28.25 


24.8 


MDM 1.3 m 


20.80±0.2 


2002 Oct 7.17 


33.8 


Palomar 200" 


20.73±0.17 


2003 Jul 2.47 


302.1 


Palomar 200" 


21.00±0.45 
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Table 2. Radio Observations of the afterglow of XRF 020903. 



Date 


At 


1.5 GHz 


4.9 GHz 


8.5 GHz 


22.5 GHz 


(UT) 


(days) 


(/^Jy) 


(/^Jy) 


(/^Jy) 


(/^Jy) 


2002 Sep 27.22 


23.8 






1058±19 




2002 Sep 29.11 


25.7 


197±70 


552±42 


746±37 


590±125 


2002 Sep 30.23^ 


26.8 






892±166 




2002 Oct 2.06 


28.6 


210±91 


788±45 


765±41 


0±270 


2002 Oct 10.13 


36.7 


294±91 


832±47 


640±40 




2002 Nov 13.02 


74.0 


175±56 


380±38 


190±30 




2003 Mar 3.71 


181.3 




215±59 






2003 Mar 7.72 


185.3 




112±45 






2003 Mar 14.74 


192.3 




127±32 






2003 Mar 17.69 


195.3 




103±40 






2003 Sep 15.22 


376.8 


43±66 


21±34 







VLBA observation 



